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RESUMEN 

Se ha incorporado al sistema de computacion CEINCI-LAB programas 
que permiten calcular la curva de capacidad sismica resistente en porticos pianos 
con aisladores sismicos de triple pendulo de friccion FPT, y, encontrar la 
respuesta sismica mediante el Metodo del Espectro de Capacidad. 

La curva de capacidad sismica que relaciona el cortante basal con el 
desplazamiento lateral maximo de la estructura se hallo aplicando la tecnica del 
pushover multimodal. El comportamiento no lineal del aislador se lo realizo con el 
modelo de tres fases presentado por Me Vitty y Constantinou (2015). 

Se presenta un manual de uso de los programas de CEINCI-LAB para 
encontrar el punto de desempeno utilizando el espectro de la Norma Ecuatoriana 
de la Construccion NEC-15. Como aplicacion se analiza el desempeno sismico de 
un portico, del Bloque Estructural 1, del nuevo Centro de Investigaciones de la 
Universidad de Fuerzas Armadas ESPE. Los resultados se comparan con los que 
se hallan al aplicar el Metodo Espectral. 

Palabras Claves: Sistema de computacion CEINCI-LAB. Tecnica del Pushover 
multimodal. Metodo del Espectro de Capacidad 



ABSTRACT 

It has been incorporated to the CEINCI-LAB computer system, programs 
that allow to calculate the resistance capacity curve in plane frames with triple 
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friction pendulum isolators FTP to find the seismic response through the Capacity 
Spectrum Method. 

The seismic capacity curve that relates the base shear with the structure 
maximum lateral displacement was found using the multimodal pushover 
technique. The nonlinear behavior of the isolator was done with the three phases 
model provided by Me Vitty and Constantinou (2015). 

It is presented a training manual for use the CEINCI - LAB programs to 
find the performance point according to the spectrum of the Ecuadorian 
Construction Norm NEC-15. As an application it is analyzed the seismic 
performance of a frame that is part of the structural block number 1 of the new 
“Fuerzas Armadas” University (ESPE) Research Center. The results are compared 
with those found using the spectrum method. 

Keywords: CEINCI-LAB computer system. Multimodal pushover technique. 
Capacity Spectrum Method. 

1. INTRODUCTION 

Una de las formas de visualizar el comportamiento sismico que tendra una 
estructura es encontrar la curva de capacidad sfsmica resistente que relaciona el 
cortante basal V, con del desplazamiento lateral maximo D t , mediante analisis 
estatico no lineal, mejor conocida como Pushover en ingles y Tecnica del Empujon 
en castellano. Chopra y Goel (1999, 2001, 2002). 

Cuando se encuentra la curva de capacidad sfsmica resistente, se halla la 
secuencia como las diferentes secciones de una estructura van ingresando al 
rango no lineal, hasta cuando se forman las rotulas plasticas. Esta visualization da 
la posibilidad de mejorar un diseno estructural si las rotulas plasticas se forman, 
por ejemplo, en el nudo inicial de las columnas del primer piso; si esto sucede el 
usuario tiene la posibilidad de incrementar las secciones o la armadura para que 
las rotulas plasticas se formen primero en las vigas de los pisos superiores. Barbat 
etal. (2015). 

Para emplear la tecnica del pushover es necesario, definir las relaciones: 
momento curvatura, momento rotation o fuerza desplazamiento, en las diferentes 
secciones de los elementos, que definen el comportamiento no lineal de los 
elementos que conduce al calculo de la rigidez a flexion, corte o caga axial. En 
este contexto el ASCE/SEI 41 presenta tablas que facilitan el calculo de los 
diagramas indicados para elementos de hormigon armado y acero. Estas tablas 
han sido incorporadas al CEINCI-LAB. Aguiar et al. (2015). 

Con la curva de capacidad sfsmica resistente, se pasa al espectro de 
capacidad de la estructura que relaciona el desplazamiento con la aceleracion 
espectral y en esta grafica se coloca el espectro de demanda sismico, en el mismo 
formato y al solo visualizar las dos curvas se infiere el desempeho sismico de la 
estructura. 

Pero con las dos graficas: espectro de capacidad de la estructura y 
espectro de demanda sismico, se determina el punto de demanda o punto de 
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desempeiio, aplicando el Metodo del Espectro de Capacidad, que fue desarrollado 
por Freeman et al. (1975), Freeman (1978). 

Existen varias formas de realizar el pushover con el sistema de 
computacion CEINCI-LAB , y son las denominadas: Monotonica, Multimodal y 
Cfclica. Mora y Aguiar (201 5), Aguiar et al. (201 6). 

El objetivo principal del artfculo es ilustrar la forma de uso del sistema de 
computacion CEINCI-LAB para hallar la curva de capacidad sfsmica resistente y el 
punto de desempeho, aplicando el Metodo del Espectro de Capacidad en una 
estructura con aisladores sfsmicos de triple pendulo de friccion FTP. 

Posteriormente se aplicara el programa para encontrar la curva de 
capacidad y el punto de demanda en un portico del Bloque Estructural 1, en dos 
condiciones, la primera suponiendo que el portico no tiene aisladores sfsmicos y la 
segunda con aisladores sfsmicos. Para los dos casos se hallara el punto de 
desempeho ante el espectro NEC-15; teniendo en cuenta que la estructura se 
halla en Sangolquf en un perfil de suelo tipo C. Bonifaz (2015). 



2. MODELO CONSTITUTIVO DEL FPT 

En la figura 1 , se presenta la nomenclatura de un aislador de triple pendulo 
de friccion, definido por doce parametros geometricos y por cuatro coeficientes de 
roce. Es bastante comun que los radios y R 4 , sean iguales; que los radios del 
aislador interior tambien lo sean R 2 = R 3 \ que la altura h 2 = h 3 , h 1 = h 4 . Que la 
capacidad de desplazamiento d ± = d A , y que d 2 = d 3 \ con esto los parametros 
geometricos se reducen a seis. 




Figura 1 Nomenclatura del aislador de triple pendulo. 
Fuente: Fenz y Constantinou (2007) 



Por otro lado, es frecuente que los coeficientes de roce jUi = ju 4 y que n 2 = n 3 . 
Con lo que se tiene 2 coeficientes diferentes de rozamiento. Para esta condition y 
para el caso de los seis parametros geometricos. Me Vitty y Constantinou (2015) 
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han definido un modelo constitutive de 3 fases, que se considera para modelar el 
comportamiento no lineal de un aislador de la tercera generation. 



En la primera fase, el deslizamiento se da en los radios de curvatura R 2 y R 3 , 
ante sismos de pequena magnitud que logran activar al aislador; en la segunda 
fase, la base deslizante rfgida llega a los topes del aislador interior y este en su 
conjunto empieza a deslizarse en las placas exteriores, como se observa en la 
figura 2. 




Desplazamiento 



Figura 2 Desempeno de aislador en Regimen II. 
Fuente: Me Vitty y Constantinou (2015) 



En el programa desarrollado, se considera el diagrama de histeresis 
indicado a la derecha de la figura 2, como que fuera bilineal. Donde el punto de 
fluencia F y esta definido por la fuerza maxima y desplazamiento maximo de la fase 
1 y el punto de maxima fuerza F u por la maxima fuerza y desplazamiento de la 
fase 2. De tal manera que las ecuaciones que definen el modelo bilineal, son: 



q y = 2 - n 2 ) R 2iBff 



F = 
y 2 R 



W 



2 ,eff 



q + jU 2 W 



F, = 



q u q y “I - 2 d.^ 
W 



2R % , 



eff 



(q - u *) + Hi W 



Hi, eff = Rt-ht Para I de 1 a 4. 



d 



I heff 

Ft 



Donde q y ,F y ,q u ,F u , son los desplazamientos y fuerzas del punto de 
fluencia y del punto que se considera de maxima capacidad del aislador. Es 
importante destacar que la fase 3 no se considera en el estudio y que corresponde 
al caso de que el aislador interior llega a los topes exteriores, con un mayor 
desplazamiento lateral; de tal manera que se esta trabajando con un modelo 
conservador debido a que el desplazamiento ultimo se considera el 
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desplazamiento maximo de la fase 2 y no el de la fase 3. Falta indicar que W es el 
peso total que Mega al aislador; q es el desplazamiento. 



Para el analisis no lineal, estatico se ha considerado dos casos, el primero 
trabajando con las rigideces elastica K e y con la rigidez plastica K p , que se indican 
en la figura 3, a este caso se denomina rigidez tangente. 



K e 



& 

Ry 



K„ 



Fu-Fy 
Ru ~ Ry 




Figura 3 Modelo constitutive bilineal para el comportamiento del aislador FPT. 
En el segundo caso se trabaja con la rigidez secante, con la rigidez K ef 



3. SISTEMA DE COMPUTACION CEINCI-LAB 

Para quienes no conocen el sistema de computacion CEINCI-LAB , se les 
recomienda la lectura del capitulo 15 del libro: Analisis Matricial de Estructuras con 
CEINCI-LAB, Aguiar (2014). Ahi veran la forma de resolver un portico o armadura 
plana y en capitulos posteriores como se encuentra la matriz de rigidez lateral y la 
matriz de rigidez en coordenadas de piso. Para profundizar un poco mas se les 
sugiere leer el libro: Analisis Dinamico de Estructuras con CEINCI-LAB, Aguiar 
(2012) y alii encontraran tambien detallado el espectro de diseno de la Norma 
Ecuatoriana de la Construccion. 

A continuacion se indica la forma de trabajar con los programas de 
CEINCI-LAB, para encontrar en primer lugar la curva de capacidad sismica 
resistente de un portico piano con y sin aisladores FPT, mediante un pushover 
multimodal y luego como se obtiene la respuesta sismica empleando el Metodo del 
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Espectro de Capacidad. Las unidades de trabajo, tienen que ser consistentes, se 
recomienda trabajar en toneladas y metros. 



3.1 Definicion de las secciones de los elementos: vigas, 
columnas y aisladores FPT 

Como primer paso se deben definir las secciones de los elementos 
estructu rales: vigas, columnas y aisladores sismicos en el caso de existir. Para la 
definicion de las secciones se debera indicar un codigo numerico para el tipo de 
material y forma de la seccion transversal, dichos codigos se indican en la tabla 1 . 



Tabla 1 Codigos de material y perfiles tipo. 





TIPO DE 
MATERIAL 


SECCION TRANSVERSAL 


Perfil 
tipo 1 


Perfil 

tipo 

H 


Perfil 

Tubular 

Rectangular 


Maciza 

rectangular 


Maciza 
tipo L 


Maciza 
tipo T 


HORMIGON 


2 


- 


- 


- 


1 


1 


1 


ACERO 


1 


1 


2 


3 


- 


- 


- 


AISLADORES 

SISMICOS 

FTP 


1.3 


1 



3.2 Vigas y columnas de acero 

Para definir una seccion de acero, sea una viga o columna, se utilizaran 
los campos indicados en la tabla 2. En la figura 2 se indica las secciones 
transversales consideradas, con su respectiva nomenclatura. 




Figura 4 Secciones transversales consideradas, para elementos de acero. 
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En el nuevo Centro de Investigaciones Cientfficas de la ESPE, las 
columnas principales son tubulares cuadradas de 550 mm., de lado y un espesor 
de 20 mm. Las vigas son tipo ‘I” de 550/300/25/12.5 (d, bs=bi, ts=ti, tw). El acero 
utilizado es A572. La forma de identificar a estas dos secciones de acuerdo al 
formato de las tablas 1 y 2, es la siguiente. 

SecC=[ 1,35150, 0,20430000, 3, .55, . 02, .55, .55, .02, .02] ; 

SecV= [1,35150, 0,20430000,1, .55, .0125, .3, .3, .025, .025] ; 

Al vector que contiene la informacion de la columna se ha denominado 
SecC, y al de una viga SecV. Se debe indicar la seccion en el nudo inicial, centra 
de luz y nudo final. 

Es importante que todos los vectores que definan la geometrfa de las 
secciones de los elementos tengan el mismo numero de columnas, si no es asi se 
colocan ceros. 

3.3 Vigas y columnas de hormigon 

La informacion para definir las secciones de las vigas y columnas de los 
elementos de hormigon armado, se indican en la tabla 3 y en la figura 5. De igual 
manera se indican tres vectores de datos por cada elemento, para el nudo inicial, 
centra de luz y nudo final. 



Tabla 3 Datos para describir una seccion de una viga o columna de hormigon 



Nombre 
de la 
seccion 


= 


Tipo de 
material 


fc 


Ey 


fy 


Tipo de 
seccion 


h 


b 


bw 


tf 


A definir 
por el 
usuario 




Segun 
Tabla 1 


Resistencia 
del hormigon 


Modulo de 
Young del 
acero de 
refuerzo 


Resistencia 
del acero de 
refuerzo 


Segun 
tabla 1 


Alto de 
la 

seccion 


Ancho 
de la 
seccion 


Espesor 
del alma 


Espesor 
de la 
alas 



d’ 


Assup 


Asjnf 


Stirrup 

separation 


Av 


Asl 


Recubrimiento o alto 
efectivo a compresion 


Area de acero 
superior 


Area de 
acero inferior 


Separation de 
estribos 


Area de las 
ramas de 
estribos 


Area de acero en el medio de la 
seccion en el caso de columnas 





As_sup 

Asl 

As2 

Asjnf 



Figura 5 Nomenclatura de vigas (izquierda) y columnas (derecha) de hormigon 

armado. 



28 



Roberto Aguiar, David Mora y Evelyn Tipanluisa 



3.4 Aisladores de triple pendulo 



Los programas de CEINCI-LAB fueron desarrollados en esta fase de 
estudio unicamente para aisladores sfsmicos tipo Triple Pendulo de Friccion FTP. 
El aislador sfsmico vendra definido por la informacion indicada en la tabla 4, donde 
se ha indicado el significado de las variables a excepcion de L^Lj que son las 
distancias desde el centra de gravedad del aislador hasta el punto medio de la 
viga inferior o superior, como se ilustra en la figura 6. 



Tabla 4 Informacion para los aisladores de triple 



pendulo de friccion. 



Nombre 
de la 
seccion 


= 


Tipo de 
material 


Tipo de 
seccion 


fy 


gy 


fu 


qu 


Li 


Lj 


kv 


A definir 
por el 
usuario 




Segun 
Tabla 1 


Segun 
tabla 1 


Limite 

elastico 


Desplazamiento 

elastico 


Limite 

plastico 


Desplazamiento 
asociado al limite 
plastico 






Rigidez 
vertical del 
aislador 




Figura 6 Modelo de un aislador como un elemento corto, de longitud L t + 



Lj 



Sea k 0 la matriz de rigidez del elemento aislador y sea k 3 la matriz de 
rigidez en coordenadas globales del elemento aislador, las mismas que se hallan 
de la siguiente manera: 



T = r-i o i t 10 // 

0 -10 0 10 . 

k 3 = T l k 0 T 

De tal manera que las distancias l u lj, sirven para hallar la matriz 
cinematica T y con ella encontrar la matriz de rigidez en coordenadas globales. El 
significado de la variables de la matriz k 0 , es el siguiente: k h es la rigidez 
horizontal; si se trabaja con las rigideces tangentes, valdra k e para el rango 
elastico y k p para el rango plastico; en cambio sf se trabaja con la rigidez secante 
k h = k ef . La rigidez vertical k v se halla con la siguiente ecuacion. 



A = 



n b\ 



E A 

k v = 7} 
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Donde: b 4 es el diametro del aislador interior del FPT; E es el modulo de 
elasticidad del acero; H es la altura total del aislador FPT. 

3.5 Asignacion de las secciones a los elementos de la 
estructura 

Una vez asignada la numeration a vigas y columnas de la estructura se 
precede a asignar las secciones antes definidas a cada uno de los elementos. 
Esta asignacion se realiza en “Tablel” de la siguiente manera: 

Tablel = 

Elemento 
1 



Elemento 

n 



Seed on en el nude 
inicial de elemento 1 


Seccion en ei centre luz de 
elemento 1 


Seccion en el nudo final de 
elemento 1 








Seccion en el nudo 
inicial de elemento n 


Seccion en el centre luz de 
elemento n 


Seccion en el nudo final de 
elemento n 



Falta indicar la geometria del portico y las cargas verticales que gravitan 
sobre el mismo; esto se lo hace en forma similar a lo indicado en el Capftulo 15 de 
Aguiar (2014), sin embargo se describe brevemente los datos que se requieren. 



3.6 Generacion de los elementos de la estructura 

Para la generacion de la matriz de rigidez y del vector de cargas 
generalizadas los datos que se deben ingresar, se resumen a continuacion. 
Primero se debe indicar el numero de pisos np y el numero de nudos restringidos 
nr. Luego para general los elementos se debe crear una matriz con la informacion 
indicada en la tabla 5. 



Tabla 5 Datos para generar los elementos de un portico piano. 



GEN 




N-= del 


Nudo 


Nude 


N D de 


1 nore-mentc- en la 


Incnemenla de la 


Increment de la 




elemento 


inicial 


final 


elementos a 
generar 


numeration 


numeracicn en el nudo 
inicial 


numeracicn en el nudo 
final 



Para la creation de las coordenadas de los nudos, la informacion 
solicitada se indica en la tabla 6. 



Tabla 6 Datos para generar las coordenadas de los nudos. 



NCDE 


= 


N c del 
nude 


Coord. X 
del nude 
inicial 


Coord X 
del nude 
final 


N c de nudes 
a generar 


In-crementa en la 
nu meracion 


Incremental de la 
longitud en X 


Incremental de la 
longitud En Y 



3.7 Asignacion de carga 

Se identifica al elemento cargado y se coloca la carga uniforme distribuida, 
en un archivo, como se indica en la tabla 7. 
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Tabla 7 Asignacion de carga vertical. 



Fm 



N n del elemenrtB cangado 



Carga dislribuida exislEniE acluando aabra el element-n 



3.8 Subrutina Pushover_3 

Es importante que antes de analizar una estructura, se revise la subrutina 
Pushover_3, y se vea las variables: 

• DeltaV, que es el incremento de cortante que se va a aplicar en cada ciclo 
de carga. Si es un valor muy bajo demandara demasiado tiempo de 
calculo y si es un valor muy alto el calculo sera muy rapido pero no es 
preciso. El incremento de cortante se reparte en cada uno de los pisos de 
acuerdo al modo de participacion y al factor de participacion modal. 

• maxderg, es otra variable a controlar, ya que es un valor que el usuario se 
impone, dependiendo del criterio de colapso adoptado. Un valor adecuado 
es 2 % para estructuras sin aisladores. Para estructuras con aisladores 
depende de la capacidad de desplazamiento de los aisladores, esta 
cantidad se divide para la altura del elemento aislador y se multiplica por 
1 00, debido a que se esta trabajando en porcentaje. 

4. ESTRUCTURA SIN AISLADORES 

4.1 Descripcion del portico 

Se analiza el portico B, del Bloque Estructural 1, del nuevo Centro de 
Investigaciones Cientfficas de la Universidad de Fuerzas Armadas ESPE, 
suponiendo que la estructura no tiene aisladores sismicos, esto unicamente para 
comparar el desempeno de la estructura sin aisladores con la estructura con 
aisladores. 

La carga muerta considerada para el analisis es D = 0.8 T/m 2 y la carga 
viva L = 0.25 T/m 2 . Para el analisis sfsmico se considera una carga de D + 0.25 L. 
En la figura 7 se indica la geometrfa del portico, con la carga vertical actuante, que 
por facilidad se considero igual en los tres pisos. En la parte inferior de la figura 7 
se indica la seccion de la columna y viga tipo; en el apartado 3.2 se presento sus 
propiedades. 
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Figura 7 Geometria del portico B, que se analiza, sin aisladores sismicos 



4.2 Curva de capacidad sismica resistente 

Se considero como criterio de colapso que la deriva global es del 2%. Con 
esta acotacion en la figura 8 se presenta la numeracion de los nudos y elementos 
con los que se realizo el calculo. 

Se utilizo un pushover multimodal y el desempeno de los elementos 
estructurales se indica a la izquierda de la figura 9; el color verde significa que la 
seccion supero el momento de fluencia y el color rojo significa que la seccion llego 
al momento ultimo. En el ejemplo ningun elemento llego al momento ultimo, 
unicamente superaron el momento de fluencia, eso significa dafio en la estructura, 
hasta llegar a un desplazamiento lateral que esta alrededor de los 35 cm. Se 
recuerda que se estan aplicando cargas laterales de acuerdo a la forma de los 
modos. 
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10 04] 11 05] 12 




Figura 8 Numeracion de nudos y elementos de portico B, sin aisladores sfsmicos 




Sistema tradicional Cortante-Desplazamiento 




Figura 9 Desempeno y curva de capacidad sismica resistente del Portico B, del 
Bloque Estructural 1 , considerando que no tiene aisladores sfsmicos. 



A la derecha de la figura 9 se presenta la curva de capacidad sismica 
resistente, considerando el criterio de colapso ya indicado; la Ifnea entrecortada 
que sobresale corresponde a la contribucion del primer modo; la Ifnea entrecortada 
que aparece en la parte inferior a la contribucion de los modos 2 y 3. La Ifnea 
continua es la resultante que se halla en cada incremento de carga al aplicar el 
criterio de combinacion modal CQC, Combinacion Cuadratica Completa (Chopra 
2001). 



Para encontrar los modos de vibracion, factores de participacion modal, 
etc., se trabaja con las coordenadas laterales, el portico es de 3 pisos, luego se 
tienen 3 grados de libertad. 
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4.3 Metodo del Espectro de Capacidad 

A la izquierda de la figura 10 se presenta la curva de capacidad sfsmica en 
el formato Desplazamiento en el tope de la estructura (tercer piso) vs el cortante 
basal y a la derecha se tiene la misma curva de capacidad, con Ifnea negra, pero 
en el formato Desplazamiento vs Aceleracion espectral. 

En la figura de la derecha se ha colocado el espectro de diseno, que 
reporta el NEC-15, para el sitio del Proyecto en un Perfil de suelo C, con Ifnea roja. 
A partir del espectro de capacidad y de demanda se calcula en forma iterativa el 
punto de desempeno aplicando el Metodo del Espectro de Capacidad. 



Sistema tradicional Aceleracion-Periodo 




Sistema de respuesta Aceleracion-Desplazamiento 




Figura 10 Curva de Capacidad Sfsmica y Metodo del Espectro de Capacidad 



En la figura 10 se observa que el punto de demanda es S d = 12.5 cm, esto 
es en un sistema de un grado de libertad, con este valor se halla que el 
desplazamiento lateral maxi mo en el tercer piso es D t = 15.6 cm, y las derivas de 
piso que se encuentran son las indicadas en la figura 1 1 . Se aprecia que la deriva 
de piso maxima es 1.31%, valor que se considera un poco alto para un Centro de 
Investigaciones Cientfficas, por los contenidos (equipos) de los mismos. 

Cuando se utiliza el Metodo Espectral para hallar la respuesta sfsmica se 
debe imponer un valor de reduction de respuesta de las fuerzas sfsmicas R, que 
es bastante incierto y la incertidumbre crece aun mas cuando se encuentran los 
desplazamientos inelasticos. Por esto es que el Metodo del Espectro de 
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Capacidad es mas exacto que el Metodo Espectral porque el metodo determina 
en base a los espectros de capacidad y de demanda el amortiguamiento 
equivalente del sistema en base al cual se reduce el espectro elastico hasta hallar 
el punto de desempeno. Adicionalmente reporta la aceleracion del punto de 
desempeno con la cual se puede hallar las perdidas en los contenidos. 



•~r r t- 



•~r T 



I I I I I I I I 

0.85 0.9 0.95 1 1.05 1.1 1.15 1.2 1.25 1.3 1.35 

Storey Drift (%) 

Figura 11 Derivas de piso halladas en el Portico B, sin aisladores 
sfsmicos. 



5. ESTRUCTURA CON AISLADORES 



En el portico B, del Bloque 1 , se van a colocar sobre los pedestales de los 
plintos el aislador de triple pendulo de friccion indicado en la figura 1 2. 




Figura 12 Geometria de aislador de triple pendulo de friccion. 

FTP8833/1 2-1 2/8-6 
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Los coeficientes de race utilizados, son: ^ = |U 4 = 0.07; |U 2 = P -3 = 0.01. 
En la estructura con aisladores sfsmicos se va a construir sobre los aisladores una 
losa de aislacion de las mismas caracterfsticas de la losa del primer piso, con esta 
acotacion, las cargas que se tienen son las indicadas en la figura 13. 




SECCION TIPICA SECCION TIPICA 

PARA COLUMNAS PARA VIGAS 



Figura 13 Cargas actuantes en la estructura con aisladores sfsmicos. 



La rigidez elastica, de pos fluencia y efectiva del aislador se calculan en 
funcion de la carga vertical que gravita sobra cada uno de los aisladores, por esta 
razon se han identificado cada uno de ellos del 1 , al 3. El 1 corresponde al aislador 
de la columna izquierda de la figura 13. Logicamente es funcion tambien de la 
geometrfa del aislador, figura 12, y de los coeficientes de race. Con toda esta 
informacion en la tabla 8 se presenta, la carga vertical total que gravita en cada 
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aislador, las fuerzas y desplazamientos que definen el modelo bilineal indicado en 
la figura 3. 



Tabla 8 Carga vertical y valores de fuerza y desplazamiento que definen el modelo 

bilineal del aislador. 



Aislador 


Carga 
vertical 
actuante (T) 


Li mite 
elastic o, Fy 

m 


Desplazamiento 
ma xi mo pa ra 
Unite elastic o, 
SUttm] 


Limite 
plastic o F Flj 
IT) 


Desplazamiento 
maxi mo para 
limite plastico., 

gu(nn) 


1 


73.095 


5.12 


0.0275 


13.34 


0.5366 


2 


136.83 


9.53 


0.0275 


25.92 


0.5366 


3 


63.79 


4.47 


0.0275 


12.03 


0.5366 



5.1 Curva de Capacidad Sismica Resistente 

Se empleo el pushover multimodal, en la estructura con aisladores 
sismicos y se trabajo con dos modelos para el comportamiento no lineal del 
aislador, en el primero de ellos se calculo con la rigidez elastica y de post fluencia; 
y en el segundo modelo con la rigidez secante. 




Figura 14 Desempeno de la estructura ante Pushover Multimodal. 
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En la figura 14 se aprecia el desempeno de la estructura ante la aplicacion 
de las cargas laterales de acuerdo a los modos de vibracion y factores de 
participation modal, Como era de esperarse unicamente los aisladores ingresan al 
rango no lineal por ese motivo estan dibujados con Ifnea verde y la estructura 
practicamente no se deforma (si existe desplazamiento lateral pero es minimo con 
respecto al desplazamiento del aislador). 

Ningun elemento de la estructura ingresa al rango no lineal, a diferencia 
del comportamiento de la estructura sin aisladores. 




Sistema tradicional Cortante-Desplazamiento 




Figura 15 Curvas de capacidad sfsmica en la parte superior se utiliza el modelo de 
rigidez tangente y en la parte inferior el modelo de rigidez secante. 

El grafico de la parte superior de la figura 15, corresponde a la curva de 
capacidad sfsmica que se obtiene con el modelo de rigidez tangente y el grafico 
que esta en la parte inferior al que se halla con el modelo de rigidez secante. 
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Como era de esperarse en el modelo de rigidez tangente es notorio el cambio de 
pendiente de la curva de capacidad, lo que no es tan evidente en el modelo de 
rigidez secante. 

Para el caso de estudio, la rigidez post fluencia es demasiado pequena por 
este motivo cuando la estructura ingresa al rango no lineal, ante un pequeno 
incremento de carga lateral el desplazamiento lateral es muy grande. Por lo que se 
recomienda trabajar con el modelo de rigidez secante. 

Si se decide trabajar con el modelo de rigidez tangente, se aconseja 
disminuir el incremento de carga lateral en cada etapa del pushover. Actualmente 
en la subrutina Pushover_3, el incremento de cortante es constante. 

En la figura 15, con linea entrecortada se presenta la contribution del 
primer modo, de los modos superiores o se observa; con linea roja la curva de 
capacidad sismica del aislador y con linea negra de la estructura completa. 

5.2 Metodo del Espectro de Capacidad 

Se recomienda mirar la subrutina: krigidez_isolator, antes de ejecutar el 
programa para ver si esta programado con modelo de rigidez tangente o con 
modelo de rigidez secante, uno de los dos modelos debe encontrarse con 
comentarios. 



Sistema tradicional Aceleracion-Periodo 




Sistema de respuesta Aceleracion-Desplazamiento 




Figura 16 Curva de capacidad sismica y Metodo del Espectro de Capacidad, en 
estructura con aisladores sismicos. 
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A la izquierda de la figura 16 se muestra la curva de capacidad sfsmica, 
hallada con el modelo de rigidez secante y a la derecha se presenta: el espectro 
de demanda sfsmico con Ifnea roja, el espectro de capacidad de la estructura con 
Ifnea negra y el punto de desempeno que tiene un desplazamiento lateral de 42.5 
cm. 



En la figura 17 se muestra la deriva de piso, en la estructura con 
aisladores sfsmicos se aprecia que existe un valor unicamente en los aisladores y 
que la deriva es nula en todos los pisos. 




Figura 17 Deriva de piso en estructura con aisladores sfsmicos. 



5.3 Uso de CEINCI-LAB 

Se presenta el programa, con el cual se halla la curva de capacidad 
sfsmica y el punto de desempeno, aplicando el metodo del espectro de capacidad, 
en la estructura con aisladores sfsmicos; se deja al lector que en base a este 
programa elabore otro para la estructura sin aisladores. 

En la figura 18 se indica la numeracion de nudos y elementos. Para los 
nudos se debe numerar primero los restringidos y para los elementos es 
conveniente que el elemento aislador se numere al final. 

Para el analisis sfsmico se considera que las vigas son axialmente rfgidas, 
con lo que se tiene un grado de libertad horizontal por piso, los mismos que se 
numeran primero, como se observa en la figura 19; luego se numeran los 
restantes grados de libertad. Se halla la matriz de rigidez de la estructura con 
todos los grados de libertad y se condensa a las coordenadas laterales. 

El grado de libertad 1, corresponde al desplazamiento horizontal del 
sistema de aislacion. 
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Figura 19 Grados de libertad considerados en el analisis. 

• Listado del programa 

1 % METODO DEL ESPECTRO CAP AC I DAD- DEMANDA 

2 % Aplicado al BLOQUE B CON AISLADORES 

3 % Dr. Roberto Aguiar, Ing. David Mora, Ing. Evelyn Tipanluiza 

4 % Nov de 2015 

5 % 

6 

7 %NOTA: TODOS LOS VECTORES DE SECCIONES DEBEN SER DE LA MISMA DIMENSION 

8 % PARA COMPLETAR LA DIMENSIONS SE PUEDE COMPLETAR CON CEROS 

9 % Unidades: toneladas y metros 

10 

11 - clear 

12- set (gcf, ' color ', 'w' ) ; 

13 - Ey=2 0430000 ; %Modulo de elasticidad del acero 

14 

15 % DEFINCION DE LAS SECCIONES 

16 

17- SecC=[l, 35150, 0,Ey, 3, .55, .02, .55, .55, .02, .02] ; %columnas rectangulares tubulares de acero 
ie - SecV=[l, 35150, 0,Ey,l, .55, .0125, .3, .3, .025, .025]; %vigas tipo I de acero 

19 - Seclsl= [1.3, 1,5. 12, 0.0275, 13. 84, 0.5366, .45, .45,4750000,0,0] ; 

20 - Secls2= [1.3, 1,9. 58, 0.0275, 25. 92, 0.5366, .45, .45,4750000,0,0] ; 

21- Secls3=[l. 3, 1,4. 47, 0.0275, 12. 08, 0.5366, .45, .45,4750000,0,0] ; %aisladores FTP 
22 

23 % AS IGNAC ION DE LAS SECCIONES A LOS ELE MENTOS DE LA ESTR.- 

24 

25 - Tablel= [ SecV SecV SecV %Eleml 

26 SecV SecV SecV %Eiem2 

27 SecC SecC SecC %Eiem3 

28 SecC SecC SecC %Elem4 

29 SecC SecC SecC £EIem5 

30 SecV SecV SecV %Elem6 

31 SecV SecV SecV %Elem7 

32 SecC SecC SecC £ElemS 

33 SecC SecC SecC %EIem9 

34 SecC SecC SecC %Eleml0 

35 SecV SecV SecV %Elemll 

36 SecV SecV SecV %Eleml2 

37 SecC SecC SecC %Eleml3 

38 SecC SecC SecC %Eleml4 

39 SecC SecC SecC £Eleml5 

40 SecV SecV SecV %Eleml6 

41 SecV SecV SecV %Eleml7 

42 Seclsl Seclsl Seclsl %Eleml3 

43 5ecls2 Secls2 5ecls2 %Eleml9 

44 Secls3 Secls3 Secls3] ; %EIem20 

45 

46 % GENE RAC ION DE LA GEOMETRIA DE LA E STRUCT . 

47 - np=4; drainer o de pisos 

48 - nr=3 ; %nume r o de mudos restringidos 
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50 


%Generacioii de los elementos de la estructura 


51 - 


GEN= [ 1 4 


5 3 5 3 3 


52 


2 5 


6 3 5 3 3 


53 


3 4 


7 2 111 


54 


3 7 


10 2 1 1 1 


55 


13 


10 13 2 1 1 1 


56 


13 


142 111]; 


57 


%Generaci6m de los rnidos de la estructura 


58 - 


NODE=[l 


0 0 1 1 11 0 


59 


3 


20.6 00000 


60 


4 


0 0.9 3 3 0 4.94 


61 


5 


11 0.9 3 3 0 4.94 


62 


6 


20.6 0.9 3 3 0 4.94] ; 


63 


% 


ASIGNACXON DE CARGAS 


64 - 


Fm= [ 1 


4.43 


65 


2 


4.43 


66 


6 


4.43 


67 


7 


4.43 


68 


11 


4.43 


69 


12 


4.43 


70 


16 


4.43 


71 


17 


4.43] ; 



6. METODO ESPECTRAL 

Con el proposito de comparar los resultados hallados con el Metodo del 
Espectro de Capacidad se realizo el analisis sfsmico empleando el Metodo de 
superposition modal, con un valor de reduction de las fuerzas sfsmicas R = 1, es 
decir se trabaja con el espectro elastico para la estructura. 

En base a la geometrfa de los aisladores, coeficiente de race y carga 
vertical actuante, se hallo que el factor de amortiguamiento equivalente % eq = 
0.2028. Con este valor se hallo el factor B por el cual se debe dividir el espectro 
para los perfodos altos en que va a trabajar el sistema de aislacion. 




De tal manera que el espectro del NEC-15, para la zona de perfodos 
cortos es el elastico ya que R = 1, y para la zona de perfodos largos esta dividido 
para B. 



Como se realizo un analisis sfsmico piano, la carga vertical se hallo en 
base a un ancho cooperante que va hasta la mitad de la luz en sentido 
perpendicular, de tal manera que esta sobre mayorada la carga vertical indicada 
en la figura 13. A mas de ello en el ultimo piso se tiene menor carga. De tal forma 
que la matriz de masas esta mayorada y la estructura es mas rfgida de lo que se 
halla en el calculo pero tanto en el Metodo del Espectro de Capacidad como en el 
Metodo Espectral se trabajo con la misma matriz de masas. 

En la tabla 9, se presentan los perfodos de vibracion, destacando que el 
perfodo fundamental es menor por la razon indicada en el parrafo anterior. Se 
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indican los modos de vibracion, el factor de participacion modal y las aceleraciones 
espectrales. 



Tabla 9 Propiedades dinamicas, factores de participacion modal y aceleracion 
espectral 





Modo 1 


Modo 2 


Modo 3 


Modo 4 


Periodo 

(s.) 


4.000 


0.594 


0.243 


0.147 


Modos de 
Vibracion 


0.156 


0.221 


-0.168 


0.078 


0.162 


0.086 


0.161 


-0.219 


0.167 


-0069 


0.165 


0.218 


0.170 


-0.216 


-0.161 


-0.077 


Factor de 
Participacion 
Modal 


6.099 


0.196 


0.027 


0.005 


Aceleracion 

Espectral 

(m/s 2 ) 


1.082 


11.091 


11.666 


1 1 .666 



Los desplazamientos laterales maximos probables que se encuentran 
aplicando el criterio de combinacion modal CQC, se indican en la tabla 10. 
Tambien se indica la fuera lateral y la aceleracion en cada uno de los pisos. 



Tabla 10 Desplazamientos laterales maximos probabl 


es 




Aislacion 


Piso 1 


Piso 2 


Piso 3 


Desplazamiento 

Lateral 

(cm) 


41.6 


43.4 


44.6 


45.6 


Fuerza Lateral 
(T.) 


11.58 


10.66 


10.82 


12.34 


Aceleracion de 
Piso 


1.13 


1.09 


1.11 


1.22 
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DERIVA POR PISO 




Figura 20 Derivas de piso encontradas en el metodo modal espectral. 

Con los desplazamientos laterales de la tabla 10, se encontro las derivas 
de piso indicadas en la figura 20, son muy parecidas a las encontradas con el 
Metodo del Espectro de Capacidad, al igual que el desplazamiento lateral. De esta 
manera se valida los nuevos programas desarrollados, sin embargo de ello a 
futuro se realizaran mas comparaciones, por lo que puede considerarse que los 
nuevos programas corresponden a la version 1 .0 

7. CONCLUSIONES 

Al sistema de computacion CEINCI-LAB se han incorporado nuevas 
subrutinas que permiten encontrar la curva de capacidad sfsmica resistente de 
porticos pianos con aisladores de triple pendulo de friccion y hallar la respuesta 
sfsmica empleando el Metodo del Espectro de Capacidad, con cualquier espectro. 
Los resultados han sido validados con el Metodo Espectral. 

Se utilizo el modelo constitutive de tres fases, para el comportamiento no 
lineal del aislador FPT, es muy sencillo trabajar con las tres fases, pero para ser 
conservadores se programo solo con las dos primeras fases, se deja la tercera 
fase como un factor de seguridad de los aisladores. 

Como ejemplo de aplicacion se encontro la curva de capacidad sfsmica 
resistente del Portico B, del Bloque Estructural 1, del nuevo Centro de 
Investigaciones Cientfficas de la Universidad de la Fuerzas Armadas ESPE y se 
aprecia que la estructura trabaja en el rango elastico, son los aisladores los que 
ingresan al rango no lineal, como era de esperarse. Se trabajo con el aislador 
FTP8833/1 2- 12/8-6 que se va a utilizar. 

Se hallo tambien el punto de desempeno y las derivas de piso, del portico 
indicado, con el Metodo del Espectro de Capacidad y con el Modal Espectral, 
encontrando una muy buena aproximacion en los resultados pero lo importante es 
ver solo el aislador es el que disipa energfa. 
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